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Experimental results for the spreading and centerline velocity decay rates for round, compressible 
jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is 
determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results 
for the far asymptotic region show that both spreading and centerline velocity decay rates, when 
nondimensionalized by parameters at the nozzle exit, decrease with increasing “jet Mach number’' 
Mj. Dimensional analysis with the assumption of momentum conservation, together with 
compressible flow calculations for the conditions at the nozzle exit, predict this Mach number 
dependence well. The analysis also demonstrates that an increase in the “potential core length” of 
the jet occurring with increasing Mj , a commonly observed trend, is largely accounted for simply 
by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing 
efficiency with increasing compressibility is inferred to contribute only partially to the latter trend. 
[SI 070-663 1(98)01 91 0-2] 


I. INTRODUCTION 

An increase in jet spreading is desired in many techno- 
logical applications, and thus, methods to achieve that are the 
subject of numerous previous investigations. Use of vortex 
generators (e.g., Refs. 1 and 2), asymmetric nozzle shapes 
(e.g., Refs. 3-5) and controlled excitation (e.g., Refs. 6 and 
7) have been among the most extensively studied techniques 
for this purpose. It should be recognized, however, that an 
increased spreading rate can be achieved only in the devel- 
oping region of a jet. Regardless of the events in the devel- 
oping region, all submerged jets should attain similar spread- 
ing characteristics far downstream. There, when the flow has 
become incompressible and the jet cross section has become 
round the velocity profiles are self-similar. Law of similarity 
dictates that a suitable quantity characterizing the jet spread- 
ing, e.g., streamwise mass flow rate, when normalized by 
local length and velocity scales, is a constant. Dimensional 
reasoning also dictates that the rate of change of the flow rate 
with streamwise distance becomes a constant. 8 

While a constant asymptotic slope for the jet spreading 
curve would be obtained under proper nondimensionaliza- 
tion, conventional way of nondimensionalization results in a 
slope that is observed to vary with jet Mach number. This 
observation, which is the focus of the paper, needs to be 
clarified first. For the problem under consideration, the prac- 
tical and most convenient choice of the parameters for data 
nondimensionalization are those occurring at the nozzle exit. 
The streamwise mass flow rate (referred to here simply as 
“mass flux”), m, is almost invariably nondimensionalized 
by the flux at the nozzle exit, m e . Similarly, the streamwise 
distance, x, is typically nondimensionalized by the nozzle 
diameter, d e . Figure 1 schematically shows commonly ob- 
served data trends under these coordinates. The curve in the 
middle is representative of a “normal,” incompressible. 


round jet. The developing region extends several diameters 
from the nozzle, as marked by the thin vertical line, after 
which the asymptotic state is reached that is characterized by 
a linear variation of the mass flux. The asymptotic slope, 
d(mfm e )/d(xfd e ), measured in many previous experiments, 
is found to be —0.28. 6,8,9 Mixing enhancement in this jet by 
a suitable technique can result in a rapid increase in the 
fluxes, typically accompanied by a lengthening of the devel- 
oping region. This is represented by the curve on the top. 
However, once the asymptotic state is reached the curve be- 
comes parallel to that of the normal jet. 

For initially compressible jets, on the other hand, first, it 
takes a longer distance to reach the asymptotic state. The 
length of the “potential core” is observed to increase with 
increasing jet Mach number, M j . (The parameter “Af } , ” 
defined clearly in the next section, represents the Mach num- 
ber at the nozzle exit had the flow expanded fully. It is used 
as the primary variable in the present study). Second, as will 
be shown clearly in the following, the asymptotic slope of 
the flux curve becomes smaller with increasing Mj . These 
trends are represented by the curve at the bottom of Fig. 1 . In 
a similar manner, the decay rate of the centerline velocity, 
U c , also exhibits a dependence on Mj . The asymptotic de- 
cay rate of U c , when expressed as 3(U e IU c )ld(xld e ) t de- 
creases with increasing M j . 

The apparent jet Mach number dependence of the 
asymptotic spreading characteristics of compressible jets is 
addressed in this paper. An understanding of these effects 
and trends is important in experimental studies dealing with 
mixing enhancement in order to have a clear perspective. It 
should also be important in turbulence modeling for com- 
pressible flows, and, perhaps, critical in numerical code vali- 
dation with jet flows. It ought to be apparent that an increas- 
ing density at the nozzle exit with increasing M } must 
influence the observed trends. As will be shown in the fol- 
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FIG. 1. Schematic of mass flux variation with streamwise distance for vari- 
ous conditions. 

lowing, only density variation does not explain the trends 
fully. An unbalanced pressure al the exit (different from am- 
bient pressure), in imperfectly expanded supersonic condi- 
tions, also comes into play. The objective of the present pa- 
per is to examine and provide an analysis for these effects. 

11. EXPERIMENTAL PROCEDURE 

The experiment was conducted in an open jet facility. 
Compressed air was supplied through one end of the plenum 
chamber fitted with flow conditioning units. The flow dis- 
charged through the nozzle, fitted to the other end of the 
plenum chamber, into the ambient of the laboratory. All jets 
involved “cold” flows, i.e., the total temperature was ap- 
proximately constant throughout the flow. Further descrip- 
tion of the facility can be found in Refs. 4 and 5. Data from 
two round nozzles are presented in this paper. One, the same 
as used for the “tab” studies reported in Ref. 2, is a conver- 
gent nozzle as shown in Fig. 2(a). The interior of this nozzle 
contracts from a diameter of 3.81 cm to the exit diameter of 
1.27 cm, the contour following a third-order polynomial. The 
nozzle has 1.27 cm long cylindrical sections at the entrance 
and the exit so that the flow enters and exits axially. 

The second nozzle, shown in Fig. 2(b), is a convergent- 
divergent one with throat and exit diameters of 2.289 and 
2.54 cm, respectively. Its interior is contoured approximately 
following the method of characteristics for a design Mach 
number, M D - 1.58. However, experiments with this nozzle 
have indicated that the flow is not completely shock free 
even at the design condition. From one dimensional nozzle 
flow analysis, based on the exit-to-throat diameter ratio, the 
following flow regimes are expected with this nozzle. (I) 
Subsonic flow throughout for M y <0.57, (II) normal shock 
existing in the diverging section for 0.57 < My <0.79, (I II) 
overexpanded flow for 0.79<Afy< 1.58, and (IV) underex- 
panded flow for 1 .58< Mj . As stated before, M j is the Mach 
number at the nozzle exit had the flow expanded fully. It is 
given by the pressure ratio through the equation Mj 
= (((/>f/p a ) (r-l)/r ~ 1)[ 2/(7~l)]) 1/2 , where p, is the total 
pressure in the plenum chamber, p a is the ambient pressure 
and y is the ratio of specific heats. Values of p,lp a corre- 
sponding to M j — 0.57, 0.79, and 1.58 are 1.24, 1.50, and 
4.11, respectively. 
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FIG. 2. Schematics of the two round nozzles, dimensions are in cm. (a) 
Convergent nozzle, (b) convergent- divergent nozzle. 

All data were obtained by Pitot probe surveys under au- 
tomated computer control. A rake of three probes was used 
to reduce data acquisition time. The probes were made of 
hypodermic needles with a 0.05 cm internal diameter and 
mounted on a streamlined “sting.” The support structure 
and the traversing mechanism were sturdy enough so that 
there was no visually perceptible vibration or yielding of the 
probes under the largest dynamic load. All data were ac- 
quired far enough downstream where the flow became fully 
subsonic in order to avoid measurement errors typical of su- 
personic flows. The distribution of Pitot pressure, on the 
cross-sectional plane at a given x, was used to calculate the 
streamwise mass flux, m = JJpUdydz, where U is the 
streamwise velocity, p is the density, and y and z are the 
transverse coordinates. In the calculation, it was assumed 
that the total temperature remained the same throughout the 
flow field and the static pressure equaled the ambient pres- 
sure. Since the integrated flux was sensitive to small mea- 
surement errors, care had to be taken in the data acquisition. 
One set of pressure transducers was employed to cover re- 
gions where the centerline Mach number was greater than 
about 0.3. A second set of transducers with higher sensitivity 
was employed in regions involving lower Mach numbers. 
Transducer zero errors were monitored before each run. 
Throughout data acquisition, the plenum pressure (p t ) and 
ambient pressure (p 0 ) were monitored and data normaliza- 
tion was done according to current conditions. Sufficient av- 
eraging time was allowed to ensure good data repeatability. 

It should be noted that there is subjectivity in the calcu- 
lation of the mass flux. Because the contribution from the 
potential flow induced by the jet is infinite in the integration, 
m = jfpUdydz, a suitable criterion needs to be followed to 
discriminate the vortical flow from the surrounding potential 
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FIG. 3. Mach number contours at xld e ~ 30 for the convergent nozzle. “Jet 
Mach number” (M y ) and contour intervals are: (a) 0.97, 0.018; (b) 137, 
0.018; (c) 1.63, 0.034; (d) 1.97, 0.057. 


flow. 6 Furthermore, Pitot probe measurements involve error 
on the outer periphery of the jet due to large turbulence and 
flow angularity. In view of these difficulties, the criterion 
followed here simply involved truncating the integration 
where the measured Mach number dropped below 1% of the 
local centerline Mach number. For further discussion of 
these considerations and limitations, the reader is referred to 
Refs. 5 and 10 in which preliminary results of this investi- 
gation were reported. Most of the mass flux data were repeat- 
able within ±4%, the repeatability being better in regions of 
higher Mach number but worse in cases involving very low 
jet Mach number. 

While the mass flux data were obtained for a limited 
number of x stations, the centerline velocity measurements 
were carried out in separate runs allowing for finer spatial 
resolution. It is worth noting that the centerline velocity U c 
is a clearly defined quantity and there is no subjectivity as 
with the mass flux data. However, measurement is not trivial 
and one must account for possible misalignment of the flow 
and geometric axes. Some further details of the measure- 
ments can be found in Ref. 10. 

III. EXPERIMENTAL RESULTS 

Examples of Mach number distributions for the 1 .27 cm 
convergent nozzle, measured thirty diameters from the exit, 
are shown in Fig. 3. The jet spreading for three supersonic 
cases is compared with that for a high subsonic case. For 
each case, ten contour levels are plotted, the interval being 
10% of the maximum contour level that approximates the 
peak Mach number in the domain. This way, the jet cross- 
sectional boundary relative to the local peak velocity is 
shown. The difference in jet spreading, depending on M ; , is 
obvious. Here, in addition to the compressibility effect, dis- 
cussed further in the following, various stages of screech 



FIG. 4. Streamwise variations of normalized mass flux for indicated M j ; 
convergent nozzle. 


cause the jet to spread differently. At Mj= 1.37 in Fig. 3(b), 
the jet undergoes the “B” (flapping) mode of screech. Under 
this condition the jet goes through the most pronounced 
spreading. At Mj= 1.63 in Fig. 3(c), the jet undergoes the 
helical “C” mode of screech and the spreading is not as 
much. At the highest value of Mj in Fig. 3(d), screech am- 
plitude is negligible, and the spreading is the least among the 
four cases. The characteristics of the screech modes have 
been addressed in many previous studies, and their effect on 
jet spreading have been observed in the early experiments of 
Glass 11 and recently explored further in Ref. 5. In the present 
paper, attention is first focused on the spreading rate of the 
jets in the far asymptotic region. 

The streamwise variations of mass flux, calculated from 
data similar to those in Fig. 3, are presented in Fig. 4. The 
data are nondimensionalized by the initial flux calculated 
from operating conditions (plenum chamber and ambient), 
with the assumption of a top-hat velocity profile at the exit of 
the nozzle. (Direct measurement of the mass flow rate by an 
orifice-meter, conducted in separate experiments, 5 showed 
good agreement with the calculated m e data.) The Pitot 
probe survey was carried out far enough downstream for 
each case so that the asymptotic region, characterized by a 
linear variation of mfm eJ was reached. In general, the mea- 
surements had to be carried out farther downstream with in- 
creasing M j . An inspection of Fig. 4 reveals that the slopes 
of the flux curves decrease with increasing M } . 

The asymptotic spreading rates were determined from 
the slopes of the mass flux curves, measured in the farthest 
downstream regions. These slopes, measured directly from 
Fig. 4, are plotted in Fig. 5(a) as a function of Mj . The solid 
line represents Eq. (4) from the analysis given in the next 
section. The trend of decreasing asymptotic spreading rate 
with increasing M j is clearly exhibited by these data. 

Dimensional analysis indicates that the spreading rate 
should scale as the square-root of the air density at the nozzle 
exit. Thus, when normalized by (p fl /pj 1/2 one might expect 
the rate to become a constant and independent of M j }' n 
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FIG. 5. Asymptotic slopes of m/m e versus xld e for the convergent nozzle; 
solid line, Eq. (4). (a) Slopes as measured, (b) normalized slopes. 


This is examined in Fig. 5(b) where the same data are replot- 
ted with such normalization. One finds that the normalized 
slope is not a constant and varies across the transonic region. 
It is this ‘ 4 anomalous’ * behavior that prompted a full inves- 
tigation with data from the C-D nozzle. Corresponding 
asymptotic spreading rate data for the C-D nozzle are shown 
in Figs. 6(a) and 6(b). Here, the flow regimes determined 
from quasi-one-dimensional flow analysis (see, e.g., Ref. 
13), as mentioned in Sec. II, arc indicated by the dotted 
vertical lines. The data in Fig. 6(b) again show that normal- 
ization by the square-root of the density ratio does not yield 
a constant value of the spreading rate. These trends, and the 
analytical predictions are discussed further in the following 
sections. 

The variations of the centerline velocity for the conver- 
gent nozzle are shown in Fig. 7(a) for different My. Note 
that the inverse of U c , normalized by the nozzle exit veloc- 
ity (U e ), is plotted. Note furthermore that U e becomes a 
constant, equal to the sonic velocity, in regimes III and IV 
(for Mj> 1 with the convergent nozzle). Clearly, the 
asymptotic slopes decrease with increasing M j . Correspond- 
ing centerline velocity variations for the C-D nozzle are 
shown in Fig. 7(b), exhibiting a more complex dependence 
on Afy . 

The asymptotic decay rates of U c > obtained directly 
from Fig. 7(a), for the convergent nozzle are shown in Fig. 



FIG. 6. Data, corresponding to those in Fig. 5, for the convergent-divergent 
nozzle. 


8(a). Corresponding data for the C-D nozzle are presented in 
Fig. 8(b). The solid curves in these figures represent the ana- 
lytical prediction described shortly. As with the spreading 
rate data in Figs. 5 and 6, normalization by ( p a lp e ) m was 
investigated for the data of Fig. 8. These results are not 
shown for brevity and it should suffice to state that the nor- 
malization yielded constant slopes only through regimes I 
and II (up to My = 1 for the convergent nozzle), and varia- 
tions occurred in regimes III and IV. The reason for these 
variations becomes clear from the analysis and discussion in 
the following sections. 

IV. ANALYSIS 

It is clear that, for the far asymptotic region {xt>d e ), the 
spreading rate, K , in mim e = K(xld e ), decreases with in- 
creasing M y . Similarly, the centerline velocity decay rate C, 
in, U e IU c — C(xfd e ) also decreases with increasing Mj. 
These trends are analyzed first. 

Dimensional analysis provides (see, e.g., Ref. 8) 

m = K x F m p\ a x, (1) 

where is a constant. The asymptotic rate of change of 
momentum, F, can be equated with forces at the nozzle exit 
as 

F = A' P 'V] + {p e -p a )A'- (2) 


I 
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FIG. 7. Stream wise variations of centerline mean velocity { U c , shown as 
l) e fU c ) for different Mj ; (a) convergent nozzle, (b) convergent-divergent 
nozzle. 


Here, A represents the nozzle cross-sectional area. The sub- 
script “a” represents ambient conditions and 4 V’ represents 
conditions at the nozzle exit. Let F= £A e p e U] y so that the 
“force coefficient,” £, is given by 

f=l + (p e -p a ) f P' U l- (3) 

Substitution in Eq. (1), with m t -A t p e U e , yields, m/m e 
= K 1 € m (A e p e U 2 e ) ]/2 p l J 2 x/(A'P'U e ). With A e =nd)lA and 
K 2 = K ] {4lir) ia , one obtains 

m/m e =K 2 Z m {p a /p e ) m {x/d e ). (4) 

In a similar manner, dimensional analysis for the asymptotic 
variation of centerline velocity yields U c — C x F m p a xn x 1 . 
One can then obtain 

U t /U c =C 2 C m (pJp,) m (xfd e ), (5) 

where Cj and C 2 are constants. The coefficient on the right 
of Eq. (4), (K= K 2 € m (p a IPe) 112 )’ as stated earlier, is ap- 
proximately equal to 0.28 for incompressible jets. Similarly, 
for incompressible flow, the coefficient on the right of Eq. 
(5), (C=C 2 £~ 1/2 (p fl /p,) 1/2 ), is —0.16 (present data; see 
also, Ref. 9). These incompressible constants, i.e., K 2 
= 0.28 and C 2 = 0.16, are assumed in the analysis. For com- 
pressible jets, one needs to calculate £ and p e in order to 
predict the variations of K and C with M j . 

As stated in Sec. II, from one-dimensional nozzle flow 
analysis, the following flow regimes are expected as the pres- 
sure ratio (P//p fl ), and hence M Jy is increased: (I) Subsonic 



FIG. 8. Asymptotic slopes of U e fU c vs xld e as a function of Mj ; solid 
line, Eq. (5). (a) Convergent nozzle, (b) convergent-divergent nozzle. 

throughout, (II) normal shock in the divergent section, (III) 
overexpanded flow, and (IV) underexpanded flow. 

A. Calculation of £ 

Through regimes I and II, p e ~p a , therefore, £= 1 [Eq. 
(3)]. Between regimes II and III (normal shock at the nozzle 
exit), there is a discontinuity in the variation of p e , and 
hence in £ (discussed further in the following). In regimes III 
and IV, M t equals the design Mach number, M D , and thus, 
p e is giv en by p e fp t =(\ +[(y- 1 )/2]Af^)~ y/(y ~ Since 
U e = M D \jyRT e , (7 being the temperature and R the ideal 
gas constant), Eq. (3) provides, £= 1 +[ l/(yMp)]( 1 
-p a fpe)- Thus, £ is known throughout regimes I-IV. 

B. Calculation of p. 

In regime I, 

M t =(((p,/p a ) <r - |Vy - D^y) . 

and in regimes ID and IV, M e = M D . In these regimes p e is 
given by, p*/p,= [l +[(y- 1 )/2]Af^]“ 1/(r_ 1} ), where p t 
= p i fRT i . For regime II (normal shock within the diverging 
section), conditions at the nozzle exit for a given p t fp Q are 
found by iteration. With standard compressible flow nota- 
tions (see, e.g., Ref. 13), the mass flow rate for choked flow, 
m* = 0532A*p l fyJY t (foot-pound-second units), is calcu- 
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FIG. 9. For the convergent-divergent nozzle, ideal variations of p e , p e , and 
£ with M j . 


lated. Assuming the normal shock at a location where the 
nozzle cross-sectional area is A, M { (Mach number upstream 
of normal shock) is found from 


1 A \ 

2 1 

2 ( 

y~ 

1 \ 1 <r + l)/(y- 

\A*} 


lr+i \ 

1 2 

-«;)] 

From, 

P a_ 

(7+1 )M\ 
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Pa 

(y-l)M? + 2 




with pt\—p t , Pa is found. This in turn provides M e from, 

/ 2 \ 1/2 

Then T e is found, with T { ~ T tl = T t2 , from, T e /T t = { 1 
+ [(y— l)/2]A/^) _l . Equations, U t — M e yJyRT e and p e 
— p a lRT e , provide U € and p e . Finally, iteration is done by 
varying A until m* = A e p e U e is satisfied. This provides 
shock location as well as conditions at the nozzle exit includ- 
ing p e . Thus, p e is found throughout regimes I -IV. Varia- 
tions of p e , £ as well as p e for the C-D nozzle are shown in 
Fig. 9. These are discussed further in Sec. V. 

C. Potential core length 

From the centerline velocity data shown in Fig. 7, a 
“virtual origin” can be defined by extrapolating the 
asymptotic straight line upstream and finding its intersection 
with the U e /U c =\ line. The distance of the virtual origin 
from the nozzle (x c ) would also provide a measure of the 
potential core length ( x p ) of the jets. The latter length de- 
notes the distance from the nozzle where the centerline ve- 
locity begins to decay. Normally, x c would be larger than 
x p , since it would take some distance, after the decay of U c 
has started, for the asymptotic linear variation to ensue. The 
length x c can be formally expressed starting with Eq. (5), by 
assigning a fixed value of unity for the left hand side. Thus, 
with D 2 = 1 /C 2 , an expression for x c is obtained as 

x c /d e = D 2 € m (p a / Pe r m ■ (6) 


As with Eqs. (4) and (5), an incompressible constant is as- 
sumed for the right hand side of Eq. (6). From the data for 
the lowest subsonic case in Fig. 7(a), this constant is as- 
sumed to be seven. That is, D ( = D 2 £ m (p a l p e )~ m ) is taken 
as ld e for incompressible flow. Note that the potential core 
length for incompressible jets is about 4 d e , but, as stated in 
the foregoing, the distance of the virtual origin is somewhat 
longer. The same calculations of £ and p e , described previ- 
ously, would then provide the variation of x c ld t with M j in 
the compressible flow regime. 

V. DISCUSSION 

The analytical curves in Figs. 5 and 6 represent Eq. (4). 
The value of the slope has been matched on the left end of 
the graphs with # = 0.28. This is the approximate value of 
the slope observed in previous experiments. 6,9 This is the 
only matching that has been done in the present analysis, the 
rest following from compressible flow equations. The agree- 
ment with the data is quite good. A similar observation can 
be made for the centerline velocity decay rates in Fig. 8. 
Good comparison is observed with the predictions from Eq. 
(5). 

Inspecting the expression for f and Eq. (4), it should be 
apparent that the value of the slope, # 2 £ 1/2 
( = K(p a lp e )~ m ), for a convergent nozzle should be 0.28 
up to Mj = 1 (£= 1 in that range and #2 = 0.28). In the limit 
of large Mj, it should reach a value of 0.367 (with y= 1.4). 
The trend is approximately followed by the data [Fig. 5(b)]. 
For a C-D nozzle, K 2 £ m should be 0.28 up to the point 
when the normal shock leaves the nozzle (demarcation be- 
tween regimes II and III), and again attain the same magni- 
tude at M d (demarcation between regimes III and IV). In the 
overexpanded regime its value should be lower than the in- 
compressible constant, since p e <p a . In the limit of large 
Mj , for underexpanded flow, its value should be lower than 
the corresponding value for a convergent nozzle, and depend 
on M d . These trends are also qualitatively in agreement 
with the data presented in Fig. 6(b). 

The discontinuity in the predicted trends between re- 
gimes II and HI is due to the assumptions in the one- 
dimensional nozzle flow analysis. There occurs a discontinu- 
ity in the total pressure variation across the normal shock 
[Figs. 2(b) and 9]. In the analysis, for regime II, it is assumed 
that the total pressure upstream of the shock equals the ple- 
num chamber pressure. Downstream of the shock, the flow 
becomes subsonic and the pressure drops to a value such that 
the static pressure at the nozzle exit equals p a . Thus, up to 
the end of regime II, p e -p a . As the shock leaves the nozzle 
exit, the total pressure at the exit goes through a discontinu- 
ous variation, attaining the value that was upstream of the 
shock. The Mach number at the nozzle exit also jumps to the 
design Mach number and remains the same throughout re- 
gimes III and IV. The static pressure at the nozzle exit, given 
by the compressible flow equations discussed in the previous 
section, is at first low in the overexpanded regime. With 
increase in plenum chamber pressure, it increases and ex- 
ceeds ambient pressure when the flow becomes underex- 
panded. This idealized variation of p e for the given C-D 
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FIG. 10. Distance of ’‘virtual origin location” vs Mj for the convergent 
nozzle; solid line, Eq (6). 


nozzle can be seen in Fig. 9. It can also be seen in Fig. 9 that 
the corresponding variations of p e and £ similarly go through 
discontinuities at the end of regime II. The discontinuities in 
£ and p e result in discontinuities in the predictions. In the 
actual flow, however, such a sharp discontinuity is not pos- 
sible. The shock would not remain normal as it approached 
the nozzle exit. There would be unsteadiness. As a result the 
sharp variations would be smoothed out. However, this issue 
concerns only a narrow range of M } between regimes II and 
III and is not relevant for the rest of the discussion. 

The reason why normalization by {p a fp e ) m does not 
yield a constant value of the slopes (i.e., constant K 2 £ m or 
1/2 ) at all A/; should also be apparent by now. A con- 
stant value of these slopes would be expected only for sub- 
sonic and fully expanded supersonic conditions. For imper- 
fectly expanded conditions the pressure term appears in the 
momentum balance equation, resulting in a value of £ that is 
different from unity. This causes the observed deviations. It 
should also be obvious that proper normalization would be 
achieved only when the term f is also taken into account in 
addition to the density ratio. Specifically, a constant value of 
the slopes would be obtained when K is normalized by 
£ Xil (pJPe) m and C is normalized by £~ 1/2 (p a /p tf ) 1/2 . In 
passing, it should be noted that the assumption, p e ~p a , 
does not strictly hold even for incompressible jets (see Refs. 
9, 14, and 15). However, the effect of the deviations in p e for 
that regime is already included in the incompressible con- 
stants (#2 — 0.28 and C , 2 = 0. 16) assumed in the analysis. 

Now, let us compare the prediction from Eq. (6) with the 
data. This is shown in Fig. 10 for the convergent nozzle. A 
few points should be emphasized before further discussion. 
First, the “potential core length,” x pi for a convergent 
nozzle is not well-defined in the supersonic regime. The ve- 
locity at the nozzle exit (U e ) is sonic. This is less than the 
velocity that would have occurred if the flow expanded fully. 
As the underexpanded supersonic flow negotiates the shock- 
expansion structure downstream, there should be an overall 
increase in the centerline velocity U c before a decrease by 
turbulent diffusion ensues. This alone would explain an in- 
creasingly longer potential core with increasing M } . The 
virtual origin location, x cy however, is a clearly defined 



FIG. 11. Distance of “virtual origin location” vs Mj for fully expanded 
jets. Symbol: Present data from the convergent-divergent nozzle; Eq. (7): 
Correlation due to Professor C. K. W. Tam based on data from Refs. 18 and 
19; Eq. (8): Effect of p e on x c , present analysis. 

quantity. It is defined the same way in the experiment and 
the analysis, and, therefore, the comparison is meaningful. 
However, x c is unlikely to represent x p in a consistent man- 
ner (say, being larger by the same ratio) throughout all flow 
regimes. Second, the analysis pertains to the asymptotic re- 
gion where self-similanty has been achieved. The domain up 
to the virtual origin, however, is characteristic of the devel- 
oping region that is affected by many parameters such as 
initial condition, nozzle shape, occurrence of screech, etc. 
Another factor that comes into play is the well-known com- 
pressibility effect that reduces turbulent mixing 
efficiency. 16,17 Reduced turbulent mixing with increasing jet 
Mach number is a factor commonly thought to result in a 
longer potential core. Finally, in continuation of the second 
point, note that there is nothing in Eq. (6) to account for the 
variation of mixing efficiency in the developing region. The 
analysis does not “know” what has happened in the devel- 
oping region that may have affected the length of the poten- 
tial core. Inspecting the derivation leading to Eq. (6), it 
should be clear that is simply the intercept of the 
asymptotic straight line, passing through the origin, with the 
1 ^ ne * The length of this intercept depends on the 
asymptotic slope that in turn depends solely on the density 
and static pressure at the nozzle exit. 

Yet, it can be seen in Fig. 10 that the data fall close to 
the prediction. This is thought to be a coincidence. Here, the 
jet is subjected to conflicting effects. A reduced mixing effi- 
ciency with increasing M } tends to elongate the potential 
core. But screech, occurring with this nozzle in the super- 
sonic regime, tends to reduce the potential core. Thus, there 
is a cancellation, and only the effect of the exit conditions 
has prevailed. 

The role played by reduced mixing efficiency with in- 
creasing compressibility is now explored further. Data for 
fully expanded flow are compared in Fig. 1 1. Only one such 
data point is available from the C-D nozzle of the present 
experiment. This is shown by the solid triangular data point. 
Additional data from the literature are represented by the 
curve marked Eq. (7). The latter warrants explanation. 

Potential core length for fully expanded jets has been 
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measured, among others, by Lau et al. 18 and Lau. 1 Correla- 
tion, x ld e = 4.2+ l.lMj, was observed to represent the data 
in Ref. 18. This correlation was modified, based on addi- 
tional data from Ref. 19, by Professor C. K. W. Tam (private 
communication, see also Ref. 20), as follows: 

x p ld e =43+\2M 2 '+f(T e IT a ), 

where f{T e IT a )~ 1.2(1 -T e /T a ) for 7,/T^l.O, and 
/( 7 ’,/ 7 ’ a ) = exp(-3.2(7’ t /r a -l))— 1.0 otherwise. In order 
to compare with the present analysis, the constant on the 
right is substituted to obtain the virtual origin location (x c ) 
as 

x c ld e = 1.0+ \2M)+f(T e IT a ). (7) 

Thus, the curve on the top of Fig. 1 1, marked Eq. (7), rep- 
resents x c inferred from the experimental results of Refs. 1 8 
and 19. The data point from the present experiment falls 
somewhat lower. This could be due to a number of factors. 
The flow with the present nozzle, as mentioned earlier, was 
not perfectly expanded and there was some residual screech. 
However, this difference is not important in the context of 
the present discussion and may be ignored. 

The prediction of x c from the present analysis, for fully 
expanded flow, follows from Eq. (6). With {=1* one can 
write 

x c /d'=D 2 (p a /p e r m , Which reduces to 

/ 1 \ 1/2 

jr c /d e =7.0^ 1 + — 2 — A^e j - (8) 

where the substitution for the incompressible constant has 
been made. Prediction from Eq. (8) is also shown in Fig. 11. 
In comparing the two curves in Fig. 11, again, it should be 
borne in mind that Eq. (7) represents correlation for x c based 
on experimental data, whereas Eq. (8) is prediction taking 
only the conditions at the nozzle exit into account. Equation 
(8) does not account for events in the developing region that 
also affect * c . It should also be emphasized that the curve 
predicted by Eq. (6) in Fig. 10 applies to a given nozzle, 
whereas Eq. (8) in Fig. 1 1 applies to fully expanded flow that 
would be obtained by different nozzles at different Mach 
numbers. 

It is apparent, from the curve representing Eq. (8) in Fig. 

1 1 , that only density variation at the nozzle exit already 
causes an increase in potential core length with increasing jet 
Mach number. There is no unbalanced pressure term here. 
Effect of other factors that influence the developing region 
(nozzle shape, screech, etc.) are also minimal. At fully ex- 
panded conditions screech should be absent, and these results 
pertain to round nozzles. The factor that is still operative 
affecting the developing region is the reduction in mixing 
efficiency with increasing Mach number. Note that while Eq. 
(8) predicts a certain lengthening of x c at a given M j , the 
data [Eq. (7)] are consistently larger. It is reasonable to infer 
that the additional increase in the potential core length, i.e., 
the difference between Eqs. (7) and (8), is mainly due to the 
compressibility effect. If Eq. (8) is expanded for moderate 
Mj , one obtains, x c /<i^7.0+0.7M^. Thus, for cold jets, 
the effect of reduced mixing efficiency with increased com- 


pressibility can be viewed as essentially an increase in the 
coefficient of the second term on the right of Eq. (7), from 
0.7 to 1.2. 

We conclude by observing, once again, that the potential 
core length increases with increasing jet Mach number not 
only due to the compressibility effect that reduces turbulent 
mixing, but also simply due to changes in pressure and den- 
sity at the nozzle exit. It should also be apparent that the 
contribution from the unbalanced pressure term can be quite 
pronounced for imperfectly expanded flows. Note the differ- 
ence in the ordinate scale between Figs. 10 and 11. The 
predicted values of x c !d e at M y = 2, for example, are about 
20 and 9 in Figs. 10 and 11, respectively. This difference is 
due to the occurrence of a high pressure at the nozzle exit 
when the flow is underexpanded. With overexpanded flow, 
on the other hand, the effect would be reversed since the 
unbalanced pressure term is negative, i.e., x c would be less 
compared to that in the fully expanded case. This can be 
easily discerned from the data presented in Fig. 7(b). 

VI. CONCLUSION 

Asymptotic spreading rates and centerline velocity decay 
rates for compressible jets are examined in this paper. When 
these rates are nondimensionalized by the scales occurring at 
the nozzle exit, following usual practice, a decrease in the 
rates is observed with increasing jet Mach number. The 
variation of static pressure and density at the nozzle exit with 
varying pressure ratio explains these trends quite well. Those 
variations also partially account for the well-known phenom- 
enon of a lengthening of the potential core with increasing 
jet Mach number. A decrease in mixing efficiency with in- 
creasing compressibility also contributes to, but is not the 
main reason for, the observed lengthening of the jet potential 
core. 
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